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Q>^ I The production cross section of heavy quarks in real and virtual photon-photon 

T-H ■ collisions has been studied. The unintegrated gluon density in the photon was ob- 

, tained using the full CCFM evolution equation for the first time. The gluon density 

\ was implemented in the Monte Carlo generator CASCADE, and cross sections for 

' heavy quark production in e'^e~ collisions were calculated and compared to LEP 

. data. Also, predictions for heavy quark cross sections in e'^e~ and 77 collisions at 

\ TESLA energies are given. 

ix. 

^ ■ 1 Introduction 

> : 
•1—1 , 

^ ■ The production of heavy flavour (6 and c quarks) in 77 collisions has been studied at 

^ LEP 11112113111001110131101 for many years and will be an important field of study 

at TESLA [TTl |T2]. It has also been studied in proton-photon collisions at HERA O 
Olia Uni im HH, and in proton-proton collisions at the TEVATRON [13 123 12111221 • 

Whereas originally the measured charm cross sections were well described by the parton 
evolution according to the standard coUinear approach, the measured beauty cross sections 
were a factor 2-4 larger than predictions. It has been shown [23 123 i2S| that an improved 
description of the heavy quark fragmentation function together with new parametrizations 
of structure functions and resumming next-to-leading logarithms in NLO QCD calculations 
result in an improved description of the B meson cross section in pp and e^e~ collisions. 
Also the application of fct-factorization and an unintegrated gluon density for the proton 
gives beauty cross sections which are in good agreement with data bom. pp collisions [23I2Z1- 
The aim of this paper is to investigate if the hh cross section also can be explained in 77 
collisions by using kt- factorization and an unintegrated gluon density for the photon. 



In a previous study |2H]; using A;j-factorization, the unintegrated gluon density in the 
photon was calculated using a simplified solution to the CCFM formalism [23 E3 lUl 1121 
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proposed by [^2] • As an alternative approach also a generalization of the GBW saturation 
model [211 ES] was used to obtain the unintegrated gluon density. In the present study we 
will use the complete CCFM evolution equations to determine the unintegrated gluon den- 
sity in the photon and calculate the heavy quark cross section with help of the CASCADE 
Monte Carlo program. 

2 The structure of the photon 

It is well known that the proton not only consists of three valence quarks, but also of 
virtual gluons that hold the quarks together. These gluons can split up in quark- ant iquark 
pairs, called sea quarks, which in turn can emit gluons, see Figure^. A photon colliding 
with a proton may therefore interact either with one of the valence quarks or with one of 
the sea quarks. 

In analogy with the proton, a photon can be seen as a flux of (virtual) quarks and gluons. 
However, the photon does not consist of any valence quarks, and must first fluctuate into 
a virtual quark-antiquark pair which can emit gluons, see Figure ^p. When a photon 
interacts with a virtual quark in another photon, the former is said to resolve the latter. 
How much of the photon that is resolved depends on the resolving power (virtuality Q^) 
of the photon. So, what a photon with a low will resolve as a quark, a photon with a 
higher might resolve as a quark and a gluon which is radiated by the quark, sharing the 
original quark momentum. A photon with a very high may also resolve the soft gluon 
as having split up in a quark-antiquark pair. This means that the more we resolve the more 
quarks we see at smaller x, since each daughter particle must have less momentum than 
the mother. The density of quarks and gluons thus depends on Q^, a phenomenon called 
scaling violation. For photons, this dependence occurs already with the anomalous 
part of the splitting, i.e. the splitting of the photon into a quark-antiquark pair. The 
density of quarks and gluons is described by the density function fi{x, iJ?) which is just the 
probability that a parton of type i (quark or gluon) is carrying a longitudinal momentum 
fraction x of the photon, at a scale fi^ = Q^. The scale is generally denoted /x^, since it 
does not always have to be the virtuality of the photon that sets the scale. The structure 
function of the photon can then be written as 

F^ix,^^')=xJ2eUn^.^^') 

i 

where Cj is the electric charge of parton i and one sums over all quark flavours. Of course, 
the photons cannot interact with the gluons, since Cgiuon = 0, and does not depend on 
the gluon structure function explicitly (in LO and NLO in the DIS factorization scheme 
[36]). However, since the quark density at small x is driven by the gluon density, it is 
possible to extract fgiuon{x^iJ?). In general, the cross section that a parton i from one 
photon will interact with a parton j from another photon and form the final state X is 

fi{Xi,fi^) ■ //(Xa, yU^) ■ ^i+j^x{Xi, X2, s) (1) 

Xi X2 

where the parton densities /7 and fj describe the probability of finding partons i and j with 
momentum fractions Xi and X2 and the partonic cross section a describes the probability 
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(a) (b) 
Figure 1: The structure of the proton (a) and the photon (b). 



that the two particles i and j will create the final state X. The partonic cross section 
is proportional to the square of a Matrix Element (ME), which in turn depends on the 
masses, virtualities, couplings etc. of the particles involved. The integral is over all xi 
and X2 allowed by the kinematics. Unlike the partonic cross section, the parton densities 
cannot be calculated from first principles. Instead, experimental fits and phenomenological 
methods must be used in order to describe the parton distributions inside the photon. 

3 /ct-factorization 

All factors in equation above only depend on the longitudinal momentum fraction x 
and the scale /i^. This means that all partons are assumed to travel in the same direction 
as the incoming particle, i.e. they have no transverse momenta. It is therefore called a 
collinear approximation. However, at large energies, where small x gluons are probed, the 
transverse momenta kt of the partons are expected to be important. Therefore, the cross 
sections are factorized jHZl EH 1^31101 into a fcrdependent partonic cross section a{x, fc^, /i^), 
where the incoming partons are treated off-mass shell, and a fcrdependent parton density 
function A{x, k^^jj?). The equivalent to equation ((H) then becomes 



The function A{x,k^, fi"^) is called the (kt-) unintegrated gluon density and is related to 
the collinear one by 



where A{x, kf, /x^) describes the probability to find a gluon with a longitudinal momentum 
fraction x and a transverse momentum kt at a scale /i^. Inverting the relation Q gives 
that an unintegrated gluon density may be obtained by differentiating an integrated gluon 
density with respect to /x^, 




dxi dx2 



dkl^dkl2Aj{xi, kf^, f/)A]{x2, k^^, f/)ai+j^x{kti, kt2, xi,X2, s). (2) 



Xi X2 




(3) 




(4) 
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Note that the unintegrated gluon density obtained in this way only depends on one scale, 
since /i^ = k^. To distinguish the different parton distributions, the notation following 
IH] is used: g{x,iJi?) is the parton distribution in the coUinear approach, while Q{x,k^) 
and A{x,k^, fj,'^) are the one- and two-scale distributions in the fcrfactorization approach. 
When using off-shell matrix elements of leading order {0{as)), some next-to leading order 
(0(a^)) processes in the collinear approach are effectively included, in addition to all 
collinear leading order processes. 

4 Heavy quark production processes 

Heavy quark production in photon collisions can be explained by three mechanisms, shown 
in Figure 121 At small ^/s, the cross section is dominated by the direct part (a-c). In this 
case, the photons couple directly to the quarks, and the process is therefore not dependent 
on the quark and gluon content in the photon. Hence, the lowest, or zeroth, order process 
(a) can be calculated in pure QED, while first order corrections, including real (b) and 
virtual (c) gluon radiation, also depend on the QCD coupling constant. These corrections 
have been calculated in jl2j and were found to contribute 30% of the cross section in the 
direct case. It has also been shown that the direct part alone (including corrections) cannot 
reproduce the measurements of b and c quark production in e~^e~ , ep and pp collisions. 

In the single resolved case (d-f), one of the photons splits up in a flux of quarks and 
gluons, where one of these gluons fuse with the other photon. The next-to-leading order 
corrections are shown in (e) and (f). Since this process depends on the quark and gluon 
distribution in the photon, it cannot be calculated using perturbative QCD, but must be 
treated phenomenologically. The direct and single resolved processes are predicted to give 
equal contributions to the cross section at y/s ~ 200 GeV (s) (M) is shown 

the double resolved case, where both photons split up and create a heavy quark pair. This 
process gives a smaller contribution than the direct and single resolved processes, since it 
depends on ^ 0.014, but it becomes increasingly important with increasing energy. 

In e"'"e~ reactions there are additional contributions to heavy quark production, which 
have been discussed in jlSI- The most important contribution comes from initial or final 
state radiation of a 'y*/Z, which decays into a heavy quark pair. The cross section of 
other processes discussed in are typically several orders of magnitude smaller at linear 
collider energies. Only the diagrams shown in Figure Eh-, d, e and g have been taken into 
account in this analysis. 

5 Higher order processes 

The processes in Figure El are all of zeroth, first or second order in the strong coupling a^, 
and the matrix elements have all been calculated using perturbative QCD. However, in high 
energy collisions it is possible to have many gluon emissions, making the zeroth, first and 
even second order calculations insufficient to decribe the data. Due to the self coupling of 
the gluons, the complexity of the processes increases rapidly with increasing order, making 
perturbative calculations very difficult. Instead one uses evolution equations to generate 
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the higher order emissions. These equations describe, under certain approximations, how 
a mother parton is spht into two daughter partons, one of which is emitted whereas the 
other continues as a propagator parton. The daughter propagator parton is then spht into 
two daughters and so on. This creates a parton ladder, as illustrated in Figure El The 
general structure of an evolution equation is 

fj = fo + J fi- Pi^j,k 

where fi and fj are the parton density functions for the mother and daughter propagator 
gluons, /o is the input gluon density function and Pi^j,k is the splitting kernel describing the 
probability that the parton i is split into two partons j and k. In the following discussion, 
the evolution starts at the bottom in Figure El and continues towards the quark box. The 
possible splittings are q qg, g qq and g gg. There are a number of different 
evolution equations, each one taking different parts of the full calculation into account and 
evolving the density functions in different variables. A short presentation of three different 
approaches follows, where the first is valid for large /i^, the second is valid for small x 
and moderate /i^, and the third is valid for both these regions. The latter is of special 
importance for this paper. 



5.1 The DGLAP evolution 

The DGLAP [111113 HEl EI evolution equation is of the form 

(iln/i^ ZTT ^ Jx ^ 
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Figure 3: Parton ladder created by gluon emissions. 



where /(x, /i^) is the density of partons carrying a longitudinal momentum fraction x 
probed at a scale /i^, and Pi^j^k{z) describes the probability that a parton i is split into 
two partons j and k with a fraction z = ^ and 1 — z oi the original parton momentum, 
respectively. The probability that a gluon splits into two gluons is given by 

Pg^U^) = ^ -2 + z{l-z) + - (6) 

where the terms and ^ are called the singular terms, since they give infinite contribu- 
tions when z ^ 1 and z ^ 0, respectively. The equation ^ thus describes the probability 
change of finding a parton of type j with momentum fraction x as we increase the scale fi^. 
In the DGLAP formalism, the gluon chain in Figure El is assumed to be strongly ordered 
in virtuality 

/x^ > \kl\ . . . > \kl\ > \kl\, (7) 

where ki is the four momentum of parton i. This means that in each splitting i —>■ i + 1, j 
one can approximate kQ = kf = . . . kf = compared to /cf+i, and since kf = mf these 
partons are considered to be massless (or on-shell). It can be shown that the ordering 
in virtuality implies that also the transverse momentum of the propagator partons are 
strongly ordered (at small z) according to 

/i' > \kl\... > |4| > \kl\ (8) 

where ku = {0, kxi, kyi,0). Hence, in each splitting one can approximate k^Q = k^^ = 
. . . fc^j = 0, which means that the DGLAP approach is a collinear approximation (which is 
obvious since there is no /cr dependence in (jS))). It can be shown, that using the DGLAP 
evolution is equivalent to resum terms of the form (a^ ln(/i^))" in the expansion of the cross 
section. Hence, the DGLAP approximation is only valid at large /i^ where these terms will 
dominate. 
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5.2 The BFKL evolution 

The BFKL I49 | IHH] evolution equation resums the terms (Q;sln(^))'^ in the expansion, 
and is thus only valid at small x. It is of the form 

di^,'g{xX')-K{klh!,'). (9) 

Here, the function K is the splitting kernel equivalent to P in (jS)). The evolution is made 
in increasing In(^), since 

Xq > . . > a;^ > (10) 

has been assumed. This implies that the emitted gluons will take a large fraction of 
the propagator momentum. However, there is no ordering in k'^ or A;^, so the collinear 
approximation can not be used, and the incoming partons of the matrix elements must be 
taken off-shell (the particles can have a virtual mass). Another important consequence of 
BFKL is that unintegrated parton densities must be used in Q, i.e. they must depend on 
kt. 




5.3 The CCFM evolution 

The CCFM jSHl IHIIl IHH IHSj evolution equation is valid both at large and small x, since 
it resums terms of both the form (asln(^))" and (a^ ln(Y3;))"- This means that at large 
X the CCFM evolution will be DGLAP-like, and at small x it will be BFKL-like. The 
CCFM evolution includes angular ordering in the initial state cascade, which means that 
the emission angles of the partons with respect to the propagator increases as one moves 
towards the quark box, 

s > a > . . . 6 > ^0, (11) 
where the maximum allowed angle S is set by the hard quark box, 

P,+p,- = T(Pel+SPe2)+4- 

This is written in the Sudakov variables, where Pq, Pq, Pei and Pe2 are the four momenta 
of the produced heavy quarks and the incoming particles, respectively (see Figure Ej), T 
and TS are the positive and negative light-cone momentum fractions of the heavy quark 
pair, and Qt is the sum of the transverse momentum vectors of the heavy quark pair. The 
momenta of the emitted gluons can be written similary, 

2 

Pi = Vi{Pel + iiPe2) + Pti, ii = (12) 

svf 

where Vi = (xj_i — Xi) is the momentum fraction of the emitted gluon, pt is the transverse 
momentum of the gluon, and s = (Pgi + Pe2)^ is, as usual, the squared center of mass 
energy. Here, we have assumed that all particles are massless. The CCFM equation is 
written as 

_2 d xA{x,klf) f d(f)P{z,kf,{q/zf) , ,^ ^ 



7 



Figure 4: The Reggeization of the gluon vertex. 



where 

Qi = . = Xi_i^/sii (13) 

is the rescaled transverse momenta of the emitted gluons and Zi = In this formahsm, 
(fTT|) becomes 

Qi > Zi_iqi_i. (14) 

When 2; — i> 1 we have > gi-i, i.e. ordering in rescaled transverse momentum, which 
means that the evolution is DGLAP-like. In the limit z ^ the angular ordering gives 
no restrictions on the rescaled transverse momentum. Also, (llOj) holds because of the 
definition of z. This means that the evolution is BFKL-like. The Sudakov form factor 
As describes the probability that there are no emissions from the starting scale /Xq to the 
maximum rescaled transverse momentum q^^^. The CCFM splitting function P is defined 
as 

P,(Z, kl {q/zf) = + ^A„.(.., kl q^l (15) 

1 Z^ Zi 

which is somewhat different than the DGLAP splitting function ©. First of all, the CCFM 
splitting function only includes the singular parts of the DGLAP splitting function. The 
other difference is that there is one additional function A^^, called the non-Sudakov form 
factor, in (fTH|l . The non-Sudakov form factor originates from the fact that, in CCFM and 
BFKL, all virtual corrections in the gluon vertex are automatically taken into account, see 
Figure |3] This is called the Reggeization of the gluon vertex. 



6 Monte Carlo simulations 

As was already mentioned, higher order processes are very difficult to calculate using 
perturbative QCD. Moreover, it is not quarks that are detected in the experiments, but 
hadrons. The transition from parton level (Figure 121 for example) to hadron level takes 
place over long distances where as is large. Hence, perturbation theory is not valid here. 
To overcome this problem, Monte Carlo (MC) generators are used. For each event, these 
programs generate all the particles, their four-momenta and all the kinematic variables 
according to certain theoretical prescriptions. The higher order gluon emissions are simu- 
lated by evolving the parton ladder according to some evolution equation (DGLAP, BFKL, 
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CCFM etc.), and the transitions to hadron level can also be implemented via hadronization 
models. 

6.1 CASCADE 

The Monte Carlo generator CASCADE jjl] is based on the CCFM evolution equation and 
thus uses unintegrated parton desities and off-shell matrix elements. For technical reasons, 
a backward evolution is used, where first the hard scattering is generated and then the 
initial state cascade (the gluon ladder) is evolved from the quark box to the incoming 
particles. In this evolution only gluon emissions are treated, that is, only the splitting 
g ^ gg is considered. CASCADE also performs the hadronization, using the Lund string 
model in JETSET/PYTHIA |S21- 

7 The unintegrated gluon density 

The evolution machinery was first tested using the one-loop (DGLAP-like) evolution for 
the proton and the photon. Then, the unintegrated gluon density for the photon was 
obtained using the full CCFM evolution. The uncertainties in the evolution are discussed 
in section 17.51 Since the distributions obtained with the CCFM evolution were compared 
with the derivative of standard gluon density functions, taken from [52], a short discussion 
of the differentiation method will be given first. Alternative methods to calculate the 
unintegrated gluon density can be found in [2H1 iM| ■ 

7.1 Numerical differentiation method 

Two numerical differentiation methods were tested: DDERIV jSH] and DFRIDR jHE], which 
are based on Rombergs principle of sequence extrapolation and Richardson's deferred ap- 
proach to the limit, respectively. The first test was to differentiate the simple function 
/(x) = — ^ with respect to x. Here, both methods gave good results, but only after ad- 
justing a parameter in DDERIV. A second test was to differentiate the GRV^ [HZl density 
function with respect to the scale /i^, then to integrate over kf and compare the results with 
the original GRV distribution. This was done because we will later use equation which 
relates the unintegrated parton density function Q{x,k^) to the integrated one gi^x,^"^). 
The adaptive Gaussian method taken from [58j was used for the integration. As seen in 
Figure the DFRIDR routine reproduced the distribution succesfully, except in the high 
X region which will be discussed later. However, with the same value of the parameter 
as used above, the DDERIV approach gave unacceptable results, and only after changing 
this parameter again the obtained distribution (which is shown in the figure) started to 
resemble the original one. Because of this poor reproduction of the input distribution, and 
its sensitivity to the choice of parameter, we decided to use DFRIDR instead of DDERIV. 



iGRV-G HO (GRV for the photon) was used. 
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Figure 6: Q{x,kf) as a function of x for different values of fi = kf. When the function 
becomes negative, the differentiation method of obtaining the unintegrated gluon density 
becomes invalid. 



The deviations from the original distribution at high x is due to the fact that the 
function becomes negative in this region. This is because of the scahng violation: as 
the scale /i^ increases, smaller distances are resolved and as a consequence we resolve more 
gluons with small x (i.e. x < 0.1) and less gluons with high x. In FigureElis shown how the 
point, where the derivative becomes negative, changes with n = kt. The smallest x-value 
for this point is obtained with kt = 2 GeV, where x ~ 0.35. Since Q{x, k^) is interpreted 
as a probability, the method becomes invalid at such high x. 



7.2 Evolution of the unintegrated gluon density 

The evolution of the unintegrated gluon density was made using the CCFM equation in a 
forward evolution procedure based on the Monte Carlo technique as described in [21] • An 
initial gluon distribution is chosen as well as a starting scale and a rescaled transverse 
momentum g<j, which is related to the emission angle (see equations (fT^ and p3|l ). for 
the first emitted gluon. The default values fJ^s = = 1-4 GeV was obtained by a fit 
to the proton structure function ^IT. The evolution then starts with qmax,o = Qsi and 
kto chosen from a Gaussian distribution around to simulate a Fermi motion of the 
partons. xq is chosen from the initial gluon distribution using the Monte Carlo method. 
The branching is then repeated until the next emitted gluon would have a higher q than 
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Qmax,o- This is repeated 10'' times for each of the 50 different maximum qmax,i chosen from 
Qmaxfi = Qs < qmax,i < qmax,50 = 1800 GeV in steps of log{qmax), thus giving a distribution 
A{x,k^,q'^) in a 50 X 50 X 50 grid of log{x), log{k1) and log{qmax)- For values between 
the gridpoints the method of hnear interpolation is used. With the unintegrated parton 
distribution available, a backward evolution scheme can be used for the simulation of the 
parton emissions. 

7.3 One-loop evolution 

To make sure that the evolution procedure is working, the evolution was first made for the 
proton with a GRV^ input distribution (using the default values /x^ = = 1-4 GeV) and 
using the one-loop approximation (HOI IHH 1121 , which means that the CCFM evolution is 
reduced to the DGLAP evolution, with the difference that the CCFM one-loop evolution 
still treats the full kinematics, i.e. it takes the transverse momentum kt into account. This 
is done by putting A^^ = 1 in the splitting function ()15|) and reducing the condition fll4p 
to ordering in transverse momentum (jH)). The obtained distribution xA{x, k^, (f ) was then 
integrated over k^, 



and compared to the original GRV distribution at different scales, see Figure [3 The small 
differences are due to the fact that the splitting function only contains the singular parts, 
and that only the splitting g ^ gg is considered. Nevertheless, the distributions are still 
in good agreement which indicates that the evolution machinery is working. 

The situation for the photon is a bit different, since the structure function, in addition 
to the hadronic component, now also consists of a pointlike component which reflects the 
splitting of the photon into a quark- ant iquark pair. Also, there are no sum rules equivalent 
to those in the proton case that constrain the quark distributions in the photon. However, 
these differences do not matter, since we only use the gluon distribution and the gluon 
splitting g — >■ gg. Figure |H1 shows the evolved distribution for the photon, using the same 
parameters as in the proton case and the GRV^ as input distribution, compared to the 
input distribution at different scales. The one-loop evolved and the GRV distributions 
show a similar behaviour, although the one-loop evolved distribution is a bit steeper at 
large scales. It is not surprising that the one-loop evolution does not give equally good 
agreement in the proton and the photon case, since the same parameters have been used 
in both one-loop evolutions, while different techniques and starting scales for the proton 
and the photon were used to obtain the original GRV distributions. 

From the above we conclude that the one-loop evolution is able to reproduce the input 
distributions reasonably well, and hence that the evolution machinery is working for the 
CCFM one-loop case. A similar test cannot be done for the full CCFM evolution, since 
the obtained distributions should be different from the ones obtained with CCFM one-loop 
evolution. However, it has already been shown [5Ej that the CCFM evolution is working 
for the proton, in the sense that it can explain data, for example bb production at the 

2 GRV 98 LO was used 
3GRV-G LO was used 




(16) 
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TEVATRON, which could not be explained with other methods. Since the only difference 
between the photon and proton evolution is the starting distribution, it can be assumed 
that the CCFM evolution will work also for the photon. 

7.4 CCFM evolution 

After making sure that the parton evolution scheme gives consistent results, the uninte- 
grated gluon density for the photon was evolved using the full CCFM evolution procedure. 
The obtained distributions integrated over are compared in Figure IHl to the GRV^ dis- 
tribution, which was used as input distribution. One of the more obvious differences is 
that the CCFM evolved distribution is much higher than the original GRV. This is due to 
differences in the definition of the integrated gluon density; one can either integrate over 
d'^kt or over dk'f. The difference is a factor tt, which explains the main part of the normal- 
ization difference. We can also note that the CCFM evolved distribution is less steep than 
the original GRV. This is a consequence of that the non-Sudakov form factor A„s — ^ 
when 2; — >■ and thereby screens the ^ factor in the splitting function. This results in a 
lower splitting probability compared to the DGLAP case (where A^^ = 1) which becomes 
visible in the small x (small z) region where the factor ^ dominates the splitting function. 

In Figure ^1 is shown the unintegrated distribution as a function of x for different 
values of kt, and as a function of kf for different values of x, compared to the derivative 
of GRV. Here, /x = g = 10 GeV. Also in Figure ITUk we see that the CCFM distribution in 
general is higher (except at very high kt) and we can also see the non-Sudakov supression 
at small x. We should also remember that the derivative method becomes invalid at high 
X because '^^[^'j^'^ ^ becomes negative. This explains the big difference between the CCFM 
evolved distribution and the GRV at kt = 2 GeV, because the derivative becomes negative 
already at x ~ 0.35 (see Figure IH)). In Figure ITIIb one can see the typical exponential 
/cj-dependence of the derivative of GRV. Since the GRV distribution is only defined for 

= fc2 > 0.3 GeV^, it is replaced with ^^^^ for k^ < k% = 0.3 GeV^ We can 

also note that the probability for the gluons to have a A;^ > /i^ = 100 GeV^ in the CCFM 
approach is significant, which would not be possible using DGLAP. This is a consequence of 
the fc^-ordering in the DGLAP evolution, where kf is always less than /x^, see (jH}, and the 
non-ordering in k'f in the CCFM evolution, where the gluons can have any kinematically 
allowed value of k^. 

7.5 Uncertainties 

In the evolution of the unintegrated gluon distributions above, the GRV density was used 
as input distribution, with fig = Qs = lA GeV given from a fit to the parton density of 
the proton [SHI- A natural question is how sensitive the results are to the selected starting 
distribution and the parameters chosen. Therefore, the evolution was made also with the 
SaS^ jM] parton distribution as input (and keeping fig = = 1-4 GeV), and also with the 
GRV parametrization as input but with fig = 1.4 GeV and qg = 1.2 GeV, fig = 1.4 GeV 

''GRV-G HO was used 
^SaS-G ID was used 
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and qs = 1.6 GeV, and /is = 1 GeV and (js = 1.4 GeV. The various distributions obtained 
are compared to the previous ones in Figure ^2 and ^| When comparing the integrated 
distributions (Figure [TT|) we see that there is only a weak dependence on the starting 
value Qs, such that the distribution becomes somewhat flatter at large g^. This can be 
understood, since a larger starting angle would mean that it takes less emissions to reach 
the maximum value (since there is angular ordering), and less emissions would mean that 
the gluon entering the matrix element (gluon n in Figure would have a larger longitudinal 
momentum fraction x. A change in /ig makes a significant difference, where a smaller 
gives a distribution which is lower at small x. One can also see that the evolution depends 
on the choice of input distribution, where the main features of the input distribution 
remains after the evolution, see also Figure El 

Also when comparing the unintegrated gluon densities in Figure El we see that the 
dependence on is small while changing fis makes a significant difference, especially in 
the region 1 < kf < 10. This is probably due to details in the evolution: In order to avoid 
soft emissions a cut-off scale fig is defined. Thus if the propagator gluon has fcj > /i^ a real 
emission is allowed, but it is not if kt < ^s- Since there is no ordering in kt in the CCFM 
evolution any propagator gluon along the evolution chain might get kt < fig after a real 
emission. In that case, the evolution of the propagator is continued until kt > fis and a 
real emission is allowed. 



8 Heavy quark cross sections 

The obtained unintegrated gluon density was used as an input to the MC generator CAS- 
CADE to calculate heavy quark cross sections in e^e~ collisions, which are compared to 
LEP data. Also, predictions for heavy quark cross sections in e~^e~ and 77 collisions at 
TESLA energies are given. 



8.1 Cross section calculation test 

The e~^e~ cross sections can be calculated using the relation 

f dQl dQl 2 2 

where is the photon flux in the electron, W^.^ is the center of mass energy in the 77 
system, Qi is the virtuality of photon i and cr^^^bb.ccx is the cross section for the subprocess 
77 — >• bb, ccX. 

The direct and single resolved cross sections for b and c quark production in e~^e~ 
collisions were calculated with CASCADE and compared with the results in [42j. Since 
the direct cross section is independent of the gluon content in the photon, and hence 
independent of evolution equations, these cross sections should agree. In heavy flavour 
production, CASCADE uses off-shell matrix elements as calculated in jnZj. Errors were 
found in (eq. B16) where a factor NcCg was missing in the direct cross section (iVc = 3 
being the color factor and Cq being the EM charge of the produced quarks). With these 
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Figure 11: The integrated CCFM evolved distributions for the photon with different combi- 
nations of fj,s, Qs and input distributions. 
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Figure 12: The unintegrated CCFM evolved distributions for the photon with different 
combinations of /ig, Qs CLnd input distributions. 
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Figure 13: The GRV and SaS distributions compared at the starting scale ji — 1.4 GeV. 
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factors taken into account, the direct cross sections from CASCADE and from leading and 
next-to-leading order calculations 112) were in good agreement, see Tabled For the single 
resolved case, the DGLAP evolution was used with the gluon density according to GRV^ 
to describe the photon structure and with cuts and parameters as in |12]. As seen, also 
these cross sections where in good agreement. The aim of this comparison is just to show 
that the basic machinery for cross section calculations in CASCADE is working. 





CHARM 






(^Dmipb) CTl-RES 


{ph) 


(GeV) 


DREES CASCADE DREES 


CASCADE 


180 


296.8 (396.3) 328 269.9 (305.8) 


217 


91.2 


191.2 (254.8) 213 93.38 (109.5) 


80 




BEAUTY 






(^Dmipb) (yi-REs 


(p6) 


(GeV) 


DREES CASCADE DREES 


CASCADE 


180 


1.280 (1.582) 1.35 1.371 (1.644) 


1.19 


91.2 


0.663 (0.823) 0.71 0.302 (0.371) 


0.27 



Table 1: Direct and single resolved charm and beauty cross sections calculated with CAS- 
CADE and compared to leading and next-to-leading order (in parenthesis) calculations in 
(DREES). Here, /i^ = 2m^, A = 0.34 CeV, nit = 4.75 GeV and = 1.6 GeV In 
CASCADE, the scale factor 0.5 was used to obtain the correct scale. 



8.2 e+e cross sections at LEP 

In Figure El are shown the obtained cross sections for b and c quark production in e~^e~ 
collisions calculated with the /cj-factorization approach, using off-shell ME and a CCFM 
unintegrated gluon density evolved from a starting distribution given by GRV, compared 
to the coUinear (DGLAP) predictions and data [H I2l El El El CHI • Here, the masses nib = 
4.5,4.75,5 GeV and rric = 1.3, 1.5, 1.7 GeV were used for the CCFM cross sections, while 
nib = 4.75 GeV and nic = 1.5 GeV were used for the cross sections calculated with DGLAP. 
This variation of the quark masses gives an estimate of the uncertainty in the calculation. 
One can see that the fc^-factorization approach describes the charm data reasonably well 
and that it gives larger cross sections compared to the collinear method. One can also 
see that the choice of mass makes a significant difference, such that a smaller mass gives a 
larger cross section. Nevertheless, the cross sections obtained with 1.3 GeV< rric < 1.7 GeV 
all lie within the experimental uncertainties. 

The cross sections for bb production, on the other hand, are not well described. However, 
in the CCFM approach, one has to determine carefully the normalization of the cross 
sections. The reason for this is that the normalization of the unintegrated gluon density 
does not have to be the same as for the standard gluon densities generated according to 
the collinear scheme, since the integration is made over different regions (i.e. the integral 

^GRV-G IHO was used, as in gl] 
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Figure 14: The charm and beauty cross sections calculated with kt-factorization compared 
to the DGLAP approach and data. 
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Figure 15: The direct and resolved cross sections in heavy quark production calculated with 
CCFM and DCLAP. 



over dkl is from to /x^ in the collinear approach, whereas in the CCFM formahsm the 
integral is over all kinematically allowed k^) and the matrix elements are different. This 
normalization factor has to be applied to the single and double resolved cross sections, but 
not to the direct cross section since it is independent of the gluon density. The contribution 
of the direct and resolved cross sections are shown in Figure ^1 As seen, the differences 
between CCFM and DGLAP are due to differences in the resolved contributions of the 
cross sections. Normalizing the CCFM cross section to the charm data at ^/s = 200 GeV 
for rric = 1.5 GeV, by applying a normalization factor n = 1.7 to the single and double 
resolved cross sections, gives a much steeper behaviour of the total cross section, as seen 
in Figure ^1 However, the CCFM predictions are still in good agreement with the charm 
data at ^/s < 200 GeV. The normalization also improves the situation for the bb cross 
section, where our calculations give af^i = 4.9 pb at y/s = 200 GeV. This can be compared 
to standard collinear NLO predictions which give cr^^ ~ 4 pb jTj and results based on j2H], 
giving af^i = 1.9 pb using a KMR gluon density, and = 2.7 pb using a GBW gluon 
density at the same energy. The CCFM cross sections are of the same order but a bit 
larger than NLO standard calculations, which has also been observed in ep collisions at 
HERA. Thus, the improvement is not sufficient fully to describe the bb data, which is still 
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Figure 16: The cross sections for heavy quark production calculated with CCFM and nor- 
malized to charm data at y/s — 200 GeV with a factor n — 1.7, compared to DGLAP 
predictions and data. 
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a factor 2-3 above predictions. 

8.3 e^e~ and 77 cross sections at TESLA 

The CCFM predictions for heavy quark production in e+e~ colhsions at TESLA energies 
are shown in Figure El Here, the normahzation factor n = 1.7 has been apphed for the 
resolved cross sections. With a charm mass of rric = 1.5 GeV the CCFM predictions at 
^ = 500 (800) GeV are 

ae+e-^e+e-ccx = 2770 (4165) pb. 
For beauty, the predictions are 

cTe+e-^e+e-bbx = 17-06 (30.72) pb 

using a beauty mass of rrif, = 4.75 GeV. 

Figure IT^ shows the charm and beauty cross sections in 77 colhsions at TESLA energies 
calculated using the CCFM evolution. The normalization factor n = 1.7 has been applied 
to the resolved contributions. At a/s = 500 (800) GeV, the predicted cross sections are 

a^.y-,ccx = 269 (359) nb 

for charm and 

a^^^kbx = 4.18 (6.38) nb 

for beauty. Also here, the masses = 1.5 GeV and rrif, = 4.75 GeV were used. Figure ITHl 
also shows the resolved contributions of the cross sections. As seen, the single and double 
resolved cross sections dominate the total cross section at these energies. The reason for 
this is that the direct contribution, in contrast to the resolved ones, decreases with energy, 
and is only about 1% of the total cross section at y/s = 100 GeV. 

8.4 Uncertainties 

The cross section for heavy quark production depends on the unintegrated gluon density 
and the partonic cross section. These in turn depend on a few parameters which are not 
fixed by theory but have to be determined from fits to experimental data, and therefore 
give rise to some uncertainties in the calculations. These uncertainties are discussed below 
and shown in Figure IT^ for e^e~ collisions. 

The CCFM evolution of the unintegrated gluon density is sensitive to the input gluon 
density, the starting scale fig and the starting value for the rescaled transverse momentum 
(jg. The default values in this analysis are fj,s = (js = 1-4 GeV, obtained from a fit to 
the parton distribution in the proton jHS], with a GRV input distribution. The gluon 
densities obtained with fig = 1 GeV and Qs = 1.4 GeV, and SaS input distribution with 
/is = gs = 1-4 GeV showed the largest differences compared to the gluon density obtained 
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Figure 17: Heavy quark cross section predictions for e'^e~ collisions at TESLA energies 
calculated with CCFM and com,pared to the DGLAP predictions. For CCFM, the quark 
masses rric = 1.3,1.5,1.7 GeV and = 4.5,4.75,5 GeV were used, where the smallest 
mass gives the largest cross section. For DGLAP, rric = 1.5 GeV and rrih = 4.75 GeV. 
Also, the normalization factor n — 1.7 was applied to the resolved cross sections in the 
CCFM calculations. 
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with the default option (see Figure [TT)). Therefore, the cross sections were calculated also 
with these gluon densities to get an estimate of the error. 

The strong coupling constant as is, despite its name, not a constant, but varies 
with the energy scale /i, according to (in a first approximation) 



(33-2-n;)ln(^)' 

Here, nj is the number of quark flavours with mass less than the energy scale /z, and A is a 
constant which has to be determined by experiments. The value of A also depends on the 
order that the process is calculated at. However, it is not clear how to compare the value of 
A used in the CCFM approach with the one used in NLO DGLAP, since the first order in 
the CCFM approach includes some second order diagrams in the DGLAP approach. The 
default value was chosen to be A = 0.2 GeV, which was also used in the evolution of the 
unintegrated gluon density, but A = 0.34 GeV (used in |l2j) was also tested. 

The energy scale /i is of the order of the typical momentum transfer, but the exact 
choice of scale should be such that the higher order contributions in the perturbative 
expansion of the cross section are minimal. The standard choice is /i^ = + , where 
m is the heavy quark mass and pt is the transverse momentum of the heavy quarks in the 
center of mass frame of the heavy quark system. To show the effect of the choice of /x, the 
scale = Am? is also used here as a comparison. 

The quark masses. Since the quarks are always confined in colorless hadrons, their 
masses cannot be measured directly. Instead, the masses must be extracted from hadron 
properties. However, this can be done in various ways depending on the exact definition of 
the quark masses, and different methods lead to slightly different results. Therefore, the 
heavy quark masses are varied such that 4.5 GeV < < 5 GeV and 1.3 GeV < nic < 
1.7 GeV. This variation of the quark masses gives the largest uncertainty, and can therefore 
be used as an estimate of the total uncertainty, as seen in Figure El for example. 

In Figure IT^ is shown the cross sections for heavy quark production in e^e~ collisions, 
calculated with different choices of A, /i^ and unintegrated gluon densities. One can see 
that a change in A from A = 0.2 GeV to A = 0.34 GeV gives the largest variation of the 
cross sections, while changing and the gluon density has a smaller effect. However, all 
these uncertainties are small compared to a change in the quark masses, see Figure [T31 



9 Conclusions 

The unintegrated gluon density for the photon was obtained with the full CCFM evolution 
for the first time. The CCFM evolved gluon density for the photon was used as input 
in the MC generator CASCADE, and cross sections for heavy quark production in e^e~ 
collisions were calculated and compared to LEP data. Also, predictions for charm and 
beauty production in e~^e~ and 77 collisions at TESLA energies are given. 

In e^e~ collisions, the obtained cross sections were somewhat larger than NLO DGLAP 
predictions, but in good agreement with the charm data from LEP. The CCFM approach 
made a slight improvement compared to the standard collinear approach for the beauty 
cross sections. However, the improvement did not fully account for the discrepancy to the 
bb data. 
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Figure 19: The cross section dependence on A, /i^ and different unintegrated gluon densities. 
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